New superconducting states in the Hubbard model 
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By acting on the excited states at half-filling in the positive-!/ Hubbard model with the ?7-pairing 
operators, we found some new superconducting states in the Hubbard model away from half-Hlling. 
Such eigenstates are different from the zero-pair and one-pair states introduced by Yang. It is shown 
that the one-pair states are not the ground-states and possess ODLRO. 



PACS numbers:75.20.Hr,7I.I0.Hf,7f.27.+a,7I.55.- 

The explanation of the high temperature superconduc- 
tivity (HTSC) in doped cuprates still remains a chal- 
lenge for theoretical physics. Nevertheless, some aspects 
of the underlying physics are almost clear: the charge 
carriers must be strongly correlated, since at half-filling, 
the cuprates are antiferromagnetic insulators, and this 
can not be explained definitely by energy -band theo- 
ries. The system should be different from the conven- 
tional Fermi liquids in the normal phase. The charge 
carriers are also paired and the pair is of the d-wave sym- 
metry As the simplest model for strongly correlated 
electrons, the Hubbard model is believed to explain the 
HTSC because it can presumably incorporate the above 
key features of the physics If this is correct, the model 
away from half-filling should display superconductivity at 
low temperatures. Especially, the ground-states should 
be superconducting. For a negative coupling, Yang first 
showed for a modified Hubbard model [Q that there ex- 
ist a number of states (even at half-filling) which possess 
off-diagonal long range order (ODLRO) Q and thus are 
superconducting (the ?7-pairing states are still supercon- 
ducting if the sign of W is changed) . The ?7-pairing oper- 
ators are just the pseudo-spin operators introduced in ref- 
erence . For the generic negative-C/ Hubbard model, it 
was also proved that for sufficiently large \ U\ or a spe- 
cial lattice configuration |^ , and filling, the ground-state 
does possess ODLRO. But, experiments indicate that the 
cuprates are Mott insulator at half-filling. Therefore, the 
Hubbard model coupling should be positive. Hence, we 
should study whether there exist superconducting states 
for the positive-?/ Hubbard model, especially away from 
half-filling. 

In this Letter, we first use the pseudo-spin operators to 
generate states away from half-filling from excited states 
at half-filling. We will show that these states are different 
from the 77-pairing states introduced by Yang but they 
also possess ODLRO. Then we will prove that the one- 
pair states also possess ODLRO, although not being the 
ground-states. We first give some preliminaries. The 
generic Hubbard model is defined by the Hamiltonian 



H = 



3 



<jk>a 



(1) 



where Cj^^c^^^ are the annihilation and creation opera- 
tors for electrons at site j with spin a and < jk > 
means j, k are nearest neighboring sites (therefore the 
lattice is bipartite in the sense of Lieb |^). For simplic- 
ity, we consider the 2 dimensional square lattice A with 
A^A lattice sites. We denote lattice vectors Gxi^y Then 
the generalized pseudospin operators can be defined as 
(here, the meaning of (—1)-' is the same as e^'^'^J where 
Q (tt,^)). 



'7, 



(2) 



where a = mex + ney,m and n are integers. is the 
total number operator of the electrons with spin a. The 
pseudo-spin operators S are related to the 77-pairing op- 
erators as 



(3) 



The total spin operator S and the total momentum are 
defined as 



P = E^(^kTCkT+4xCki) = PT+Pi , 

k 



(4) 



where c,vr — N, 



-1/2 



The momentum com- 



mutes with the spin and pseudospin. Since the Hamil- 
tonian (|l]) commutes with the set (P, S, S) the generic 
model therefore enjoys SU(2)ig)U(l)g)U(l) symmetry and 
the eigenstates of the Hamiltonian can be designated by 
\E, P; s, Sz', s, Sz >, where E, P; s, Sz', s, Sz are the coore- 
sponding eigenvalues of the operators H, P, S'^,Sz, S^, Sz- 
In particular \E, P; s, Sz] s, > denote the states at half- 
filling, and the s must be an integer. According to Lieb 
Q, s = for the ground-state and it can also be shown 
that s = 1^ Q is the unique ground-state. The ex- 
cited states can be denoted by P; s, s^; s, > with 
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s 7^ 0, s 7^ 0. The partial particle-hole transformation T 
is defined by 

TctT-i = (-ircl^, fc,^f-^ = c,i (5) 
and implements the transformations TST^^ = S, and 
The Hamiltonian is thus transformed into 

fH{U)f-^ = H{-U) - U{S, + 5.) + ic/iVA (6) 

The state P; s, s^; s, > is transformed into the 
state \E,'P - 7V|Q;s,0;s,s^ >. The state \E,V - 
iVf Q; s, 0; s, Sz > accommodates Na — 2sz electrons, and 
the total spin, now equal to s satisfies s < ^N\—Sz. Since 
the z-component of the spin satisfies —^Na < Sz < ^Na, 
we arrive at the following lemma. 

Lemma: The total pseudo-spin s of the half-filling 
states \E,P; s, Sz', s,0 > for the model (1) satisfies the 
constraint s < Na- 

Since [77^, H] = —Urj'l for a = 0, while 

j 

for a ^ 0, according to Yang B, the states >:= 
ir]^)"~^ril\0 > and |F„ >:= (jrj^lO > are 2n-eiectron 
eigenstates for all a. Our first result can thus be stated 
as 

Theorem 1: Suppose P; s, s, > is an ex- 
cited state at half-filling, then the states >:= 
77"|S, P; s, s,; s, > and |$+ >:= (77+)"!^, P; s, s.; S, >, 
for n < min(s, ■^), possess ODLRO and are different 
from \Ynj^_^ ^ > and \Ynj,_^^ > for s 7^ 

Proof. First let us consider the state |<i>,^ >. It con- 
tains Na — 2n electrons and we must show that they do 
not vanish. Using the angular-momentum commutation 
relations we have 

~ S'^z - Sz = S-S+, (7) 

Therefore, 

<$ri*r >=<$ois'-^'-^.i*o >=s(5+i) (8) 

Similarly, we have 

< $-1$- > = < $;_,|S2 - - Sz\^-_^ > 

= [S{S + 1) - n{n - 1)] < <f>-_i|$,-_i >, (9) 

Therefore, if |<E>o > is normalized, we have (similarly for 

m > 



(10) 



1=1 



Therefore, |<i>^ >^ 0. Since 



(11) 



1$^ > are also eigenstates. Now let us investigate the 



ODLRO in >, i.e. limi 



i-j\^co 



For this purpose, we investigate the quantity < |S — 



f3'{^-l)Nl + ^NA 



5(5+1) 



(12) 



where we used that n = as = PNa- According to a 

theorem by Yang and Zhang |10|, there exists ODLRO 

in \^I3N A > so long as a < 1. Obviously, |$o >= 

\E;s,Sz;^Na,0 > also possess ODLRO. 

Next we show that > are different from \Yna _ > 

2 ^ 

and \Yn, > for s 7^ Since the eigenvalue of 

' -j^— n,a ' 2 ^ 

corresponding to > is s{s + 1) while the one corre- 



sponding to |y. 



> IS 



Na(Na 



1), the eigenvectors 



must be mutually orthogonal. Using 



tl - 



we have 
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— —n.a 



(^t)^-« 



J. ^ A 

= (r;t)^-' 
2 ^ 2 



'iiv^V + Vz 



-vis') 



|0> 

+ '71S2)|0> 



> 



(13) 



>. This com- 



> 



Hence > is different from \Yn. 
pletes our proof of theorem 1. 

We know that for positive J7, the energy of \Yna 

is lower than that of \Yna > hy U: Ena — 

Ena = U ( but in the thermodynamic limit, the 

energy-densities of these two states are equivalent). In 
|l|, Yang proved that \Yna > is not the ground-state 

by showing that the energy of |i^iVA_ > is equal to 

2 ^ 

the energy expectation of the (non-eigen)state >:= 

(7y^)^~"|0 > where rjc is the Cooper pair, rjc '■— 
Our second result concerns | Yna ^ ^ > and 

can be stated as 

Theorem 2: The state IFn. > is not a ground state 

for Na — 2n electrons and possesses ODLRO. 

Proof. We first show that it is not a ground-state. For 
this purpose, we consider the state 



62 >:= 



(14) 



and the energy expectation value in it. From 
[U J^j nnnji^Vc] = U'nl;we have [U J^j ^jinji, (77c)"] = 
mU{ri]j)"^, therefore 
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^wCu = U{— -n-l) 

< 02102 > 2 



(15) 



Now consider the kinetic energy: T — "^ji^ fjtjk<^j„Ckcr- 
In k-space, (m = ^^ — n — 1) 

102 >= E e'PV,^,4,...cLk„acLTCQ-PiCpTlO> 

p,ki...km 

Obviously, the k's in the sum X)ki k2 k differ from 
each other. Now consider the ki-pair. For a given p G 
BZ, the ki-pair contributes when ki ^ p,ki ^ p — Q. 
For each such ki corresponds k'j^ :— Q — ki; the k'j^-pair 
also makes contributions when k'^ ^ p, k'^ ^ p — Q, thus, 
both ki" and the k'j^-pair contribute simultaneously when 



ki ^p,±(p-Q),2Q-p 



The energy of the ki-pair and the corresponding k'j^-pair 
cancel each other because of i(k) = — i(Q — k). Thus 
only when 



ki = P, or ki = Q p 



(17) 



the ki-pair will contribute to the kinetic energy. Since 
the other k's are different from ki, the corresponding 
pairs contribute nothing to the kinetic energy, and it is 
the same as 



<02|r|02 > 
< 02 102 > 



^ < ipi\T\ipi> < ip2\T\if2 > . 



< ^2W2 > 



where >= |cL(Q_p)|c|;j_p^ctj_p^c|,-f > and |^2 > 
Therefore 

< 02|T|02 > ^y-[^(p _ Q) + ^(Q _ p) + ^(Q _ p) 



< 02 02 > 



P 



132 = 2(iVA - 3)/3i+ < Q\Tiam^v'ml\^ > (21) 

Using eq. (p^), the second term vanishes. Therefore 
/32 = 2Na{Na — 2)(iVA — 3). In general, one can obtain 
the relation 



I3n = n{NA -n- l)/3„_i 



Thus we have 



NA-l{NA-n-2)\ 



(22) 



(23) 



Therefore, we should have n < — 2 for jl^n+i^a >• Now 

let us check the behaviour of C„ :=< Yn^alCr^CriCl^Cl^lYnM > 

as \r — s\ ^ oo. We define 



(16) A. < 0\vaCr^Criv''-\v^rvClcl.ii^ 



> 



We can obtain the recurrence relations 

En = niNA + 1 - n)En-i 

Dn = n{NA -n- l)Ai-i + E^ 

Cn = n{NA -n- 3)C„-i + Ai 

Since v'^Cl^Cl-^'qllO >= we have i?„ = 0, hence 

Dn = niNA -n- l)Ai-i (24) 

By direct calculation, we know Di = (A^a— 4)e'*-i'(^'-~^=\ 
therefore 



Dn = n\{NA - 4) , 3)! »Q-(R.-R.) 

^ '{Na~ 71-2)1 

Accordingly, we can obtain 



(25) 



<P) + tip) + t{-p) + m - p) + tip)] = (18) 



iNA-5)\ 
^"-"•(iVA-4-n)!^^ 

7VA-3-n, 



where we have used that i(k) = t(— k),t(k) = — i(Q — k). 
Finally, 



< 02|g|02 > ,.,Na 

; — U n — 1) 

< 02 02 > ^ 2 ^ 



+ Dn[l + i 
+ ( 



A^A - 1 - n^ 
iVA-3-n A^A-2-n, 



(19) 



i.e., 102 > has the same energy as that of |^na_„^ >, 
but 102 > is not an eigenstate, hence I^na _„ ^ ^'^^ 
the ground-state. 

The next is to show that it possesses ODLRO. For 
this purpose, we have to know the norm of {Yn^^ ^ ^ >. 

Define 



Therefore 



A^A - 1 - n Na-u ' 
Na-S-h NA-2-n NA-'i 
Na-I-u ' Na-ti '"Na-2 



)] (26) 



(20) 



For n = 1, direct calculation gives /3i = NAiNA — 2). For 
P2, we can use the relation [77,77^] — 2rjz to obtain 



, ¥ Cn^n\j^^^^-^C^ + in-l)Dn (27) 

|r-s|^oo (AiA— 4 — n)! 

Since Ci = (A^a - 4)e^*5 (^'-^»), we have 
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Finally, we arrive at 



k-sHoo <Yn,a\yn,a> Na{Na - 2) ^ ' 

This completes our proof of Theorem 2. 

We now pass to the conclusions. First, when s ^ 0, at 
least one site is double-occupied in the half-filling state 
\E, P; s,Sz;s,0 > because r]\E, P; s, Sz;s,0 >^ 0. There- 
fore, all excited states at half-filling have at least one site 
double-occupied. Second, we have obtained the explicit 
expression of liin|r_s|^o2 < rjrrjl >. In fact, to show the 
existence of ODLRO, we can also use a theorem from 
Indeed, since 

<r„|r„> + (30) 

we also have 

<^"-°'^^-^-'^"-">=-^+nA.A^n^ (31) 
< Y Y ■> 2 

which is of order 0{N'^) when n is of oder A^a. Finally, we 
do not know whether the states |$„ > are ground-states 
or not. We do not know either which one of the states 
|$„ > and \Ynj^ ^ > has a lower energy. Though we 

do not know whether the ground-states away from half- 
filling are superconducting or not, we can make sure that 
there is a number of states which are superconducting 
away from half-filling in the positive-C/ Hubbard model. 
Of course, whether these superconducting states are re- 
sponsible for the HTSC remains unclear. 
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